Kondo effect in the presence of spin-orbit coupling 
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We study the T — Kondo physics of a spin- 1/2 impurity in a non-centrosymmetric metal with 
spin-orbit interaction. Within a simple variational approach we compute ground state properties of 
the system for an arbitrary form of spin-orbit coupling consistent with the crystal symmetry. This 
coupling produces an unscreened impurity magnetic moment and can lead to a significant change 
of the Kondo energy. We discuss implications of this finding both for dilute impurities and for 
heavy-fermion materials without inversion symmetry. 
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Introduction.- Kondo effect, i.e. screening of the 
impurity magnetic moment by the Fermi sea of itin- 
erant electrons, is one of the best-knovifn examples 
of correlations-driven phenomena in condensed matter 
physicsi. A system involving a periodic array of such 
impurities interacting with conduction electrons (the so- 
called Kondo lattice) is believed to provide a minimal 
model for heavy-fermion compounds^. Historically the 
Kondo screening was detected via resistivity measure- 
ments in dilute magnetic alloys, but recent advances in 
scanning tunneling spectroscopy allowed observation of 
this phenomenon on the atomic scale^^— and manipula- 
tion of individual Kondo resonances^. 

Details of the band structure of the host metal usu- 
ally do not qualitatively influence the Kondo ground 
state, although they affect characteristic energy scales 
of the problem such as the Kondo temperature, Tk, be- 
low which the impurity spin is screened. Similarly, in 
the presence of spin-orbit scattering when spin is not a 
good quantum number, classification of the states by par- 
ity still allows mapping of the impurity problem onto a 
Kondo model with essentially same parameters but with- 
out the local spin-orbit interaction (SOI)^, in agreement 
with experimental^. 

In non-centrosymmetric materials a distinct non-local 
(dependent on the gradients of the crystal potential) SOI 
appears. This interaction is odd in electron momentum 
and couples it to the electron spin^Q. The influence of 
this type of SOI on manifestations of the Kondo effect 
was discussed only recentlyiir— in quasi two-dimensional 
(2D) systems for specific cases of Rashba or Dresselhaus 
SOI, and in the context of topological insulator a^^d^ . 

In Ref. 11 it was concluded that, to lowest order, the 
Rashba SOI only leads to a rescaling of the electron band- 
width and leaves the Kondo temperature essentially un- 
changed. A similar verdict was reached in Ref. [l^ in 
the framework of the Anderson model for a half-filled 
/-band. However these results rely heavily on the spe- 
cific form of the Rashba SOI term and 2D single-particle 
density of states. This particular combination allows re- 
duction of the Kondo Hamiltonian with SOI (equivalent 
to a multichannel problem, see below) to a single-channel 
model without spin-orbit coupling. What happens with 



Kondo screening in more realistic and interesting cases, 
e.g. three-dimensional materials without inversion sym- 
metry or systems with a non-Rashba SOI that do not 
allow the above simplification has not been explored. 

In this Communication we consider a single spin-1/2 
impurity interacting with a system of electrons in a non- 
centrosymmetric metal at zero temperature. Due to the 
explicit inversion symmetry breaking, the single-particle 
Hamiltonian that describes the conduction band contains 
an odd in momentum spin-orbit term compatible with 
the crystal symmetryii. We determine the ground state 
properties of the resulting Kondo Hamiltonian by gener- 
alizing Yosida's variational method^^ to take into account 
the spin-orbit splitting of the Fermi surface (FS), as well 
as all values of the total spin of the electrons and the 
impurity. In contrast with previous works^'^^'^'* our anal- 
ysis is valid for any form of SOI and the electron band 
structure, and incorporates the essentially multichannel 
nature of the problem. We give general expressions for 
the Kondo binding energy and show that the SOI may 
lead to an enhancement of the Kondo effect compared 
to that of a centrosymmetric material with the same pa- 
rameters. Because the SOI breaks SU (2) symmetry, the 
impurity spin no longer forms a singlet with the Fermi sea 
and is only partially screened. This conclusion is qual- 
itatively similar to the situation in 2D helical metals^-. 
Since our goal is to investigate only effects associated 
with SOI we ignore possible spin anisotropy terms anal- 
ogous to those appearing in the study of impurities near 
sample surfacesi^,. 

We first set up the variational framework, and then 
present results for the Kondo binding energy, total spin 
in the ground state, and the impurity spin susceptibility. 

Variational formalism. - The Kondo model describes 
a localized magnetic impurity interacting with a single 
band of conduction electrons 

k 

This Hamiltonian is defined on a lattice with N sites; S 
is the impurity spin [S — 1/2) located at site iq, 
are Pauli matrices, creates a fermion at site i with 
spin a = (t,i) (4„ = Vl/^Si is its momen- 
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turn space counterpart), and Sapik) is the single-electron 
dispersion. We take Jk > 0, assume summation over re- 
peated indices, and set h= 1. 

For a single band with SOI the matrix £a^(fc) can be 
writteiiii as: 

The scalar ek is the dispersion without SOI. The latter 
enters through the real pseudovector = — r_fc which 
is determined by the point group symmetry of the crystal. 
It is convenient to diagonalize Ea/s explicitly by introduc- 
ing the helicity basis Cka — {Uk)a\dkx with A ~ ±1 and 
unitary matrix Uk such that TiJJj^crUk = cr^|rfc|. In 
this representation the band energy is diagonal, Skx — 
{Ul)\aeafj{k){Uk) px = efc + A|rfc|. Note that Tk breaks 
parity but preserves time-reversal, hence Skx — £-k.x be- 
cause of the Kramers theorem. Now we can rewrite the 
Kondo Hamiltonian as: 



H ^^Ekxdl^^dkx- 

k.X 



(1) 
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A'X'dkX. 



To understand the influence of SOI on the Kondo 
screening, we use the Yosida-like^S trial wavefunction 



AkMxOiskx - eF)\M)dl^\FS), (2) 



k,M,X 



where AkMX are variational amplitudes, M = (f, i) la- 
bels impurity states and |FS) is the filled Fermi sea 



|FS) 



n 4+ n 4jo>. 



Efc- <£f 



The Heaviside function 9{ekx — sp) limits summation to 
the energies above the Fermi level. The expectation value 
of the Hamiltonian Eq. ^ in the state of Eq. ^ is 

{i{j\H\^) = ^Al,]^,pyAkMx^i£k'X' - £F)Q{£kx- sf)^- 

X £kx5x'x5k'kSM'M + -^^M'M{Ul,(TUk) , 

with the implicit summation over all indices in the r.h.s. 
In this expression we omitted the A-independent ground 
state energy of the Fermi sea, Eq = (FS|i?|FS). Com- 
puting the expectation value of the Kondo interaction 
requires decoupling of the product 

{FS\dp,x'dl^dkpdl^\FS) = 9{ep,x, - £f)x 

X [Sal3Skk'Spp'S\\>6{eF — eka) + 

+ 5kp5k'p'5l3x5aX'0{£pX — £f)] • 

In this equation the first term has the form 
Hk^J^^^F - £ko.){UlfTUk)^^ = 2(FS|5e|FS), where 



Se = (1/2) X^/e ""a/JCfcaCfc/? IS the total clcctron spin; 
(FS|5'e|FS) ~ due to the time-reversal symmetry. 

Minimizing {ip\H\'ip) w.r.t. AkMX^ one obtains an 
eigenvalue equation: 

(sfcA — E)AkMxQ{£kx — £f) = —JKSMM'd{ekx — £f)x 
-^^0(ep^ -ei.)(C/tcrC/p)^^^v4pM'r„ (3) 



where again the summation over doubly repeated indices 
is assumed. To proceed further, we introduce 



BkMX = ^'^{Uk) ^jO{ekrj — £F)AkMr], 



'X,, 



which allows us to rewrite Eq. ([3]) in the form 
BkMa — —JrSmr^ 



(4) 
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This object plays the role of the ground-state wavefunc- 
tion for the system. Due to the 6'-function in the defini- 
tion of BkMa all fc-summations are over the entire Bril- 
louin zone. 

We shall now use Eqs. ([2]), ([3]), (01) to compute the 
Kondo energy, total spin of the system and impurity 
magnetic susceptibility in the most general form. Then 
we apply obtained expressions to several instructive ex- 
amples: (i) quasi-2D systems with symmetry C'^y (with 
Rashba or Drcsselhaus SOI), and (ii) cubic crystals with 
symmetry T ov O. 

Kondo energy.- The energy eigenvalue E in Eq. ([3]) 
is obtained by summing Eq. (|4]) over k 



k.X 



{Uk)^, 



£kX 



X0 



with X]^Ja = i^/N)J2k-^kMa- The A-dependent terms 
between two [/-matrices can be decomposed as 



di£kX 
where 



£f)5x'x/ {£kx ~ E) = Sx'xK+ + c 



A'A' 



2 V £k+ — E Ek- — E 



(5) 



Because Ukcr^Ul = Tk(T/\Tk\ and Tk is odd, while 
£kxik) is even in fc, the term containing (t^k_ does not 
contribute to the sum and we find 



X 



Ma ~ -^SMRCaP ^ K+X 



N 



(6) 



Clearly, the lowest-energy solution has the "singlet" 
structure in the helicity space: Xmo — {5M\5a- — 
i5Af^<^Q+)/v^- Then the sum is computed as 

e-E ^ 2 SE' 
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where W and ep ^ W are the half-bandwidth and Fermi 
energy respectively, SE = ep — E <^ W, and gp is the 
density of states (DOS) in the A-branch at the Fermi 
level. From this expression we finally obtain the energy 
of the Kendo bound state 



5E = We~'^^^'^'^^^i^'^p\ 



(7) 



When the SOI is absent g+ = g~ = g^p and Eq. 
reduces to the well-known result for the usual Kondo 
effecti: 5E'^°^ = M^exp[-2/3J;f5^°^]. If the characteris- 
tic SOI energy for electrons near the FS is A50 ^ Ef, we 
expand the DOS at the Fermi level up to the second or- 
der ~ gf+Ksodg^FldKso+^od^9Fl'^dKlo, where 
the derivatives are evaluated at A50 = 0. We estimate 
in a metal dgp/dKso • 



>^9f''/£f, and d'^g^/dk%o 



ifff V^F with the sign depending on the curvature of the 
DOS around the FS. Therefore {g% + gp)/2 - gf^{l ± 
^so/^f)- While in typical materials Aso/^f ^ 0.1 and 
the above correction is only ^1%, the exponential form 
of the Kondo energy, Eq. ([7]), makes the effect non- 
negligible 

SE/SE^°'> = e^^so/^lJKQp^ ^ ^±aIo/sfJi< ^ (-g) 

Assuming Aso ^ Jk this gives ~ 10% change in the 
Kondo energy relative to its value SE'^^'> without the SOI. 

It is instructive to apply the general expressions (O 
and ([8]) to two examples with parabolic bands (with an ef- 
fective mass m) and a linear in k SOI: (i) quasi-2D tetrag- 
onal systems characterized by Rashba (Dresselhaus) spin- 
orbit couphng with = ^so^ x = ^sok) 
and tetragonal axis pointing in the z direction, and (ii) 
non-centrosymmetric cubic crystals^^ with Ffc — Agofc. 
The coupling constant A50, which has units of veloc- 
ity, introduces a natural energy scale £50 — ^^%ol'^i 
and is related to Aso via A50 ~ Aso/kp where kp is 
the helicity-averaged Fermi momentum. Consequently, 

(-so ~ ^so/^F- 

In case (i) the DOS per helicity A and for positive en- 
ergies is given by gx{e > 0) = g'^p^ [I - \^/esol{(^SO + 
with g^' = m/27r. Consequently gp + gp = 'ig^p^ and 
Eq. ([7]) yields no correction to the Kondo energyiiii^: 
6E — SE^'^K This conclusion is specific solely to 2D sys- 
tems with parabolic bands and linear SOI. Of course, 
cubic in momentum SOI terms will introduce corrections 
of the form ([5]). In contrast, for case (ii) we have 



without SOI the Kondo Hamiltonian ([T]) can be reduced 
to a one-dimensional form which simply reflects the fact 
that only electrons with zero orbital angular momentum 
couple to the impurity2£. When the SOI is taken into 
account, such reduction is not always possible because of 
the ?7fc-matrices in Eq. ([T]) which entangle different or- 
bital harmonics. While in systems with Rashba SOI one 
can still decouple orbital channels by introducing suit- 
able linear combinations of c-operators and show that 
only one of them enters the Kondo terroii, other forms 
of SOI, e.g. case (ii) considered above, do not allow such 
simplification. Thus the validity of Eq. ([T]) is only re- 
stricted by the variational Ansatz ([2]). 

Total spin in the ground state.- In the standard 
Kondo problem^ at zero temperature the impurity is fully 
screened by the Fermi sea and the net spin of the sys- 
tem vanishes. This is not the case in the presence of a 
SOI. Because of the latter, even without the impurity 
the electron system has a non-zero spin, (FSjS'glFS) — 
(1/4) ^ \{Ul(TUk)^.\^eiep - SkM^k. ~ ef). This ex- 
pression is finite due to the mismatch between Fermi sur- 
faces for different helicities. Therefore, our goal in this 
part is to compute the difference between net spins in the 
Kondo and normal metal phases, (V'|(5'-|-5e)^|V')/(V'IV') — 
(FS|5'g |FS). Note that due to time-reversal symmetry of 
the problem, the total spin polarization along any direc- 
tion still vanishes. 

Using Eq. ([5]) and the singlet structure of Xua [see 
discussion after Eq. ([6|)], we can rewrite Eq. ^ as 

so that the norm of the state ([2|) becomes 

kMa kMa 

= (^)'E(4 + -^-)- 

The cross-terms ^ k^k- vanish due to the same argu- 
ment as that used in deriving Eq. ([5]). Next, we consider 
the expectation value of SSe 

{tp\SSe\tp) — - ''^ BIj^.j,^, Sm' MCTa'aBkMa = 



4 V 2 



5A(e > 0) 



m^Aso / 1 + e/2eso 



When eso ^ £f this DOS leads to an enhancement of 
the Kondo energy 6E/6E^°'> = Q'^so/epJKg^p^ ^ agree- 
ment with Eq. 0. 

It is important to emphasize that Eqs. (O and ([5]) 
correspond to a generally infinite channel Kondo prob- 
lem even in the parabolic band approximation. Indeed, 



where again there are no cross-terms and T^^'' — 

{S'S'S^)wR{a'Uk<y''Ula^Uk(7^Ula^)Yj. Since 5 is a 
spin-1/2 operator, we can evaluate T using the relations 
^*R'i''^R'R'^y-y^Ri = -^ab and crVV^ = ieuj + (SuSsj + 
SjiSis — Sis5ij)a^ where euj is the fully antisymmetric ten- 
sor: 



R'y _ 



1 



Rf 



{26_ 



R'RO-y'-y 



7'7 
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Collecting the above expressions we have 



\SS,, 



-3/4 Efc 4 + 1/4 Efe^'- 



and 



4E.(4 + «'-) ' 

with the second term in the r.h.s. coming from {ip\S'^\'ip) . 
In the absence of spin-orbit band splitting k_ = 
and the above expression implies complete screening. 
In the presence of SOI the change in the total spin is 
also finite and for cases (i) and (ii) considered above 
{^\{S + Se)"^]^) - (FS|52|FS) - A|o. In principle this 
change can be determined from local magnetic measure- 
ments, but more precise methods than the one used here 
may be needed to determine spatial dependence of the 
spin-spin correlations. 

Impurity spin susceptibility. - Finally, we consider the 
linear susceptibility of the system. Since our focus is on 
the effect of SOI, we shall make a simplifying assump- 
tion that the system is either cubic or tetragonal with 
magnetic field pointing along the c-axis, and that the 
dominant effect of the field is on the impurity spin. In 
both cases the Hamiltonian, Eq. ([I}, acquires a pertur- 
bation SH = —iJ,BSz, where fi = gfis, fJ-B is the Bohr 
magneton and g is the appropriate Lande factor. 

In order to account for SH we need to change E 
E + hM in Eqs. (g]) and © with h = and M = ±1. 

A solution is sought in the form: Xmu = ^^sijv/a + ^t^Afa 
with y (F*) the normalized singlet (triplet with zero 
total spin ^-projection) basis states. 



1 - 3P Q \ fxs 
-3Q l + P [xt 



0, 



where P = {Jk /4:N)J2k,M + 't-^'^) ^nd Q = 
(JK/4iV)X;fc,M^^'«+(-S' + '^M)- To lowest order in 
^B/5E, the ground state energy becomes E — ep^SE — 
jjp'B^ /^5E. Therefore changes in the Kondo energy ([7]) 
are straightforwardly reflected in the spin susceptibility 



X 



-d^E/2dB'^ = n^/8SE. 



Discussion.- Stimulated by the interest in non- 
centrosymmetric /-electron materials^'^'*, we investi- 
gated the influence of the lack of inversion symmetry 
on interaction between conduction and localized elec- 
trons by studying a single impurity Kondo model with a 
SOI in the conduction band. Using a simple variational 
framework-iiiSi we presented results for the ground-state 
properties of the system, valid for any form of SOI and 
band structure of the host metal, even in cases when 
one cannot reduce the problem to a single- channel Kondo 
Hamiltonian. It is the variational nature of our approach, 
what allows us to deal with a multichannel model. In par- 
ticular, we demonstrated that: (1) the SOI can lead to 
an exponential change of the Kondo temperature; (2) as 
the SOI explicitly breaks SU{2) symmetry the Fermi sea 
does not completely screen the impurity spin, allowing an 
extra magnetic degree of freedom in the Kondo phase. 

Although a similar exponential enhancement of the 
Kondo temperature was found in Ref. [l^ we note that 
their result is physically different from ours. The reason 
for this distinction is the fact that in Ref. [H the authors 
started from an Anderson model and used a Schrieffer- 
Wolff transformation^^. Although this is the usual way 
to "freeze" charge fluctuations at the impurity, in the 
presence of SOI it can lead to unexpected results, such 
as the Dzyaloshinky-Moriya coupling between impurity 
and conduction electrons spins, which appears because 
of virtual transitions of localized electrons into the con- 
duction band where they accumulate a phase due to SOI. 
On the contrary we started with a Kondo model that in- 
cludes only spin fluctuations. Thus modiflcations to the 
Kondo energy, Eq. ([T]), compared to its value in a cen- 
trosymmetric material originates purely from SOI. 

Our findings lead to an intriguing question regarding 
the infiuence of SOI on the physics of the spin-1/2 Kondo 
lattice model. It is known^'.^- that the heavy-fermion 
(Kondo screened) state competes with magnetic phases. 
In the presence of a SOI impurity spins are not com- 
pletely screenedS^ and may order, thus leading to a co- 
existence of the heavy fermion state and magnetism. We 
leave investigation of this problem for a future work. 
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